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Abstract 

We calculate the full one-loop electroweak radiative corrections to the cross 
section of single VF-boson inclusive hadroproduction at finite transverse momentum 
(pt)- This includes the 0{a) corrections to W+j production, the 0{as) corrections 
to + 7 production, and the tree-level contribution from W + j photoproduction 
with one direct or resolved photon in the initial state. We present the integrated 
cross section as a function of a minimum-pT cut as well as the pj- distribution for the 
experimental conditions at the Fermilab Tevatron and the CERN LHC and estimate 
the theoretical uncertainties. 
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1 Introduction 



The hadroproduction of single W bosons via the Drell-Yan process in pp colhsions at 
the CERN SppS led to the discovery of this particle in 1983 [T]. Nowadays, this process 
serves as a standard candle to calibrate and monitor the luminosity of hadronic collisions, 
since its cross section is rather sizeable and W bosons are straightforward to identify 
experimentally thanks to their simple and distinct decay signature. The quahty of the 
luminosity determination is thus limited by the precision to which this cross section is 
predicted theoretically. It is, therefore, mandatory to calculate higher-order radiative 
corrections. At present, they are known at next-to-leading order (NLO) [2] and next-to- 
next-to-leading order (NNLO) [3J in quantum chromodynamics (QCD) as well as at NLO 
[3] in the electroweak sector of the standard model (SM) of elementary particle physics. 

In order for the W boson to acquire finite transverse momentum (pt), it must be 
produced in association with one or more other particles. To lowest order (LO) in QCD, 
the additional particle is a gluon (g), quark (g), or antiquark (g), materialising as a hadron 
jet (j). The corresponding partonic subprocesses are of 0{aas)- Their cross sections 
are presently known at NLO [5] and NNLO [6] in QCD, i.e. at C(aa^) and 0{aa^), 
respectively. Very recently, also the one-loop electroweak corrections, at 0{a'^as), were 
considered [7|. This is also the topic of the present paper. However, as explained below, we 
actually study somewhat different cross section observables and arrange for our results to 
be manifestly infrared (IR) safe by avoiding kinematic cuts that destroy the inclusiveness 
of massless quanta. In fact, the situation is complicated by the circumstance that both 
photons and gluons can appear as bremsstrahlung. 

The observable we thus wish to investigate is the differential cross section for the in- 
clusive hadroproduction of single W bosons with finite pt at 0{a'^as). Specifically, we 
concentrate on the pt distribution and the integrated cross section as a function of a 
minimum-pr cut, leaving the distributions in other observables, such as rapidity, for fu- 
ture work. At LO, the system X recoiling against the W boson is purely hadronic, while 
at O(a^as) it can include a photon (7). We thus also have to consider W + j produc- 
tion, whose LO partonic cross sections are of C(a^), because its NLO QCD correction 
contributes at the very order we are aiming at. In fact, the real radiative corrections 
to W + j and W + 'y production receive contributions from a common set of 2 ^ 3 
partonic subprocesses. In the former (latter) case, IR singularities from the radiation of 
soft photons (gluons) cancel against similar contributions from virtual photons (gluons) 
by the Bloch-Nordsieck theorem [S]. Similarly, the IR singularities from coUinear final- 
state radiation (FSR) cancel against similar contributions from the virtual corrections by 
the Kinoshita-Lee-Nauenberg theorem [9J. The residual IR singularities from coUinear 
initial-state radiation (ISR) are factorised and absorbed at 0{a) {0{as)) into the parton 
distribution functions (PDFs). This procedure leads to manifestly IR-safe cross section 
observables and translates into unique and simple event selection criteria on the experi- 
mental side. 

By contrast, the notion of electroweak NLO corrections to W+j production comprises 
a conceptual problem. In fact, in the treatment of final-state coUinear singularities caused 
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by the parallel emission of a photon from an outgoing (anti) quark line, one is led to intro- 
duce a cut in an appropriate separation variable. Within the coUinear phase space region 
thus defined, the (anti) quark-photon system is effectively treated as one particle whose 
momentum is identified with that of j and thus subject to an acceptance cut in transverse 
momentum, pr(j) > Pt^O)' ensure the experimental observation of j. This includes 
phase space configurations where the photon essentially carries all the momentum, while 
the (anti)quark can, in principle, be arbitrarily soft. This will not generate any soft IR 
singularities. However, since (anti)quark and gluon jets can, in general, not be distin- 
guished experimentally on an event-by-event basis, the same recombination procedure 
needs to be applied to a gluon-photon system in the final state as well. This time, a soft 
gluon will inevitably produce an IR singularity, which can only be canceled by the NLO 
QCD corrections toW + 'j production, so that one falls back to the symmetric procedure 
outlined in the preceding paragraph. Formally, this soft-gluon singularity can be avoided 
by applying the cut just to the transverse momentum of the gluon, even if it is 

accompanied by a collinear photon. However, such a prescription is purely academic and 
quite unsuitable for experimental implementation because (anti)quark and gluon jets are 
treated on different footings. 

By crossing external lines, the LO partonic subprocesses of + 7 hadroproduction 
can be converted to those of W + j photoproduction with one incoming photon partici- 
pating directly in the hard scattering (direct photoproduction). The emission of photons 
off the proton can happen either elastically or inelastically, i.e. the proton stays intact 
or is destroyed, respectively. In both cases, an appropriate PDF can be evaluated in the 
Weizsacker- Williams approximation [T0llllfl2j . Since they are of 0{a), these direct pho- 
toproduction contributions are of 0{a^). Incoming photons can participate in the hard 
scattering also through their quark and gluon content, leading to resolved photoproduc- 
tion. The contributions from direct and resolved photoproduction are formally of the same 
order in the perturbative expansion. This may be understood by observing that the PDFs 
of the photon have a leading behaviour proportional to aln[M'^ / Aqq-^)) oc a/as, where M 
is the factorisation scale and Aqcd is the asymptotic scale parameter of QCD. Although 
photoproduction contributions are parametrically suppressed by a factor of a /as relative 
to the 0{a'^as) corrections discussed above, we shall include them in our analysis because 
they turn out to be quite sizeable in an extensive region of phase space. 

This paper is organised as follows. In Section [21 we list the partonic cross sections 
at LO and explain how to evaluate the hadronic cross section from them. In Section [3l 
we discuss in detail the structure of the NLO corrections. In Section HJ we present our 
numerical results for pp collisions with centre-of-mass (cm.) energy = 1.96 TeV at 
the Fermilab Tevatron and pp collisions with = 14 TeV at the CERN Large Hadron 
CoUider (LHC). 
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2 Conventions and LO results 



We consider the hadronic process 

A{pa) + B{pb)^W{p) + X, (2.1) 

where the four-momentum assignments are indicated in parentheses. We work in the 
colhnear parton model of QCD [13] with Uf = 5 massless quark flavours q = u,d,s,c,b, 
neglect the masses of the incoming hadrons, A and B, and impose the acceptance cut 
Pt > p't^ on the transverse momentum px of the W boson. (We assign masses to the 
partons 'y,g,q,q only to regulate soft and collinear IR singularities in intermediate steps 
of our calculation.) 

Specifically, denoting ui = u, U2 = c, di = d, d2 = s, and d^ = b, the relevant partonic 
subpro cesses include 

(2.2) 
(2.3) 
(2.4) 



Ui + dj - 




' + 


g, 


Ui + g- 




' + 


dj, 


dj + g- 









at 0{aas 



Ui + dj - 




' + 


7, 


tij + 7 - 




' + 


dj. 


dj+-i- 




' + 


Ui, 



(2.5) 
(2.6) 
(2.7) 



at C(a^), and 



Ui + dj^W^ + g + 7, (2.8) 
Ui + g^W^ + dj + 7, (2.9) 
dj+ g-^W^ + Ui + -i, (2.10) 

at 0{a^as)- The partonic subprocesses involving a W~ boson emerge through charge con- 
jugation. Processes fl2.2l) - fl2.5p must be treated also at one loop, Oi^a^as). Processes (12.61) 
and (12.71) contribute to direct photoproduction and processes f l2.2p - (l2.4l) to resolved pho- 
toproduction. Since photon emission off protons happens at 0{a), it is sufficient to deal 
with photoproduction at tree level. In summary, we calculate the cross section of pro- 
cess ( 12. ip at NLO as the sum 

^AB^WX ^ ^ ^ W,^^ ^ ^^.^^ ^ W,, ^ ^W,^ (2.11) 

where a^^ and (y^^^) processes (I2.2p - (l2.4p at tree level and one loop, a^^ and 

"^Gtls) "^^^ process (12.50 at tree level and one loop, a^-^"' is due to processes ( 12. 8p - 
(I2.10p at tree level, and is due to processes ( 12. 6p and (12.70 via direct photoproduction 
and due to processes (I2.2p - (l2.4p via resolved photoproduction, both at tree-level. 
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The minimum-pr cut is necessary to stay away from the regions of phase space that 
are sensitive to the colhnear IR singularities due to the q 1 /g + <f , Q 1 / 9 + T i 
1 Iq ~^ 1 '^Ti ci^nd 7/(7 — s> g + g* sphttings, which are present aheady at LO. Here, an 
asterisk marks a virtual parton. The cross section (j^^^^^^ of the hadronic process fl2.1l) 
is related to the cross sections a'^^^'^^W of the partonic subprocesses, 

a{pa) + h{p,) ^ W{p) + c{p,){+d{pd)), (2.12) 

where a, b,c,d = 7, g, q, q and pa = XaPA, Pb = XbPs with scaling parameters Xa, Xb, as the 
incoherent sum 

^AB^wx^g^p^ > p-*) = Yl f'dTCi^{T)a^'-^<'\s,pT > pTI (2.13) 

a,b,c{,d) "^^0 

where 5* = {pa+PbY and s = (pa+Pb)'^ = tS are the hadronic and partonic cm. energies, 
respectively, r = XaXb, and 

Ci^'ir) = f ^/.m(x„ M')fb,B f-, M^) (2.14) 

Jt Xa \Xa / 

is the parton luminosity defined in terms of the PDFs fa/A{xa, M"^), fb/sixb, M^). Here, 
M denotes the factorisation mass scale. Introducing the short-hand notation w = M^, 
we have 




(2.15) 



In order to obtain ^-"^'-^^^^(d)^ have to evaluate the transition matrix elements 
q-ab^w^d) q£ processes (12.121) . square them, average them over the initial-state spins and 
colours, and sum them over the final-state ones, which leads to |7"«''^w'c(d)|2 ^^le order 
of our calculation, 7"'^*^^^=^ jg calculated at tree level, while T"-'^^^'^ may receive also 
one-loop contributions, T"''^^^ = Tg"''^^^ + T^^^^m^^, so that 

^q-ab^Wf = ^j-ab^Wc^^ + 2Re [(r^c^b^Wcy j-ab^Wc-j _ (2.I6) 

Then we have to integrate over the partonic phase spaces imposing the minimum-p^ cut. 
In the following two subsections, we describe how this can be conveniently done for the 
two- and three-particle final states, respectively. 

Since we are dealing with charged-current interactions of quarks, Ui and dj, the 
Cabibbo-Kobayashi-Maskawa quark mixing matrix Vij appears. At tree level, the cross 
sections of processes ( ]2.2p -( !2rT0l) contain the overall factor and a part of the one- 

loop corrections is proportional to V*jVij'V^j,Vi'j, where Ui' and dj' are virtual quarks. 
Since we neglect all down-quark masses, we can sum over the indices of the virtual and 
outgoing down quarks to trigger the unitarity relation Yl'j=i — ^n'- ^^^e of 
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incoming down quarks, we can absorb the residual appearances of |Vijp into a redefinition 
of their PDFs, as [7J 

3 

U/a{x,M') = J2 m'f,^/A{x,M'), (2.17) 
i=i 

and similarly for down antiquarks. Therefore, it is sufficient to calculate the partonic 
cross sections for the flavour-diagonal case, with Vij = 6ij. 

2.1 Two-particle final state 

If parton d is absent in process ( I2.12p . we supplement s by two more Mandelstam variables, 
t = [pa ~pY ^"^^ — {Pb — p)"^- Four- momentum conservation implies that s + t + u = w, 
and we have = tu/s. The partonic cross section entering Eq. fl2.13p is evaluated as 

a'^'-'^%s,PT > P't^') = / dpT—^ , (2.18) 

jpcut apT 



where p™^^ = (s — w)/{2^/s) and 

^^ab^Wc 



Pt 



IT^^^^^Y + {t^u). (2.19) 



dpT Ms^{s - wY - Asp^ 

For the reader's convenience, we list the differential cross sections of processes (12.2 
(12. 7p . in the conventional form 



da 



ab^Wc 



dt IQtis' 



\'Jab^Wc\^^ (2.20) 



at LO. The Feynman diagrams contributing to processes (12.21) and (12. 5p are displayed in 
Figs. [1] (a) and (b), respectively. We have 

^^ud^w+g ^ 27raas s'^ + - 2tu 
dt ~~9si sHu ' 



, 1 + , (2.21) 

dt 12a, \ s-wj dt ^ ^ 

where = sin^^ is the sine of the weak-mixing angle. Since the VT-boson mass sets the 
renormalisation scale of the couplings, it is natural to adopt the definition of Sommerfeld's 
fine-structure constant a in terms of Fermi's constant Gp, 

a/2 

a = —Gpslw. (2.22) 

71 

The implementation of this renormalisation scheme at one loop is explained in Sec- 
tion [SUl The cross sections of processes (12. 3p . (12. 4p . (12. 6p . and (12. 7p may be obtained 
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from Eq. (I2.2ip by exploiting crossing symmetries, as 





3 


' ^da'"^^^^3' 




' dt =' 




' dt 




S<->M 










" At = 


s - 


dt 
















-3 


' 2d(3-"'^^^+T" 




' dt =■ 


^ dt 
















s - 


dt 





(2.23) 



2.1.1 Three-particle final states 



If parton d is present in process f l2.12p . then the partonic cross section entering Eq. fl2.13l) 
may be obtained through a four-fold phase-space integration along the lines of Ref . [2] . 
We work in the partonic cm. frame and choose our coordinate system so that points 
along the z direction and pd lies in the x-y plane. We denote the polar angle of by 
and the azimuthal angle of Pc by (p. As the first three independent variables, we select 
■t?, and (p, which take the values 



0<p'd< 



w 



< -i? < vr, < V9 < 27r. 



(2.24) 



In the case of process (12.81) . which contains two massless gauge bosons in the final state, 
it is convenient to take the fourth variable to be with values 



2v'- — 



1 
<^°<2 



w 



2p 



We then have 



^4:^ud^W+gj 



dp°dp°d cos -^dif 



K27r) 



^q-ud^W+g^^^g (^p^_pcnt^ _ 



:2.25) 



^2.26) 



On the other hand, in the case of processes (12. 9p and (I2.10p . which only contain one 
massless gauge boson in the final state, it is more useful to choose the fourth variable to 
be the angle ip enclosed between p^ and pd, with values 



< V < TT. 



(2.27) 



We then have 

f^'i^ug^W+d-y 



dp°d cos -^dipdip 



Pt>Pj 



P'd [y^ (v^ - W - M 
16(27r)4 [yi - 2p0 sin2(V;/2)] ' 



^j-ug^W+d^\^0 ^p^_pcnt^^ 



(2.28) 
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and similarly for process fl2.10p . In order to implement the minimum-pr cut, px needs to 
be expressed in terms of the integration variables, which is conveniently done with the 
help of Eqs. (5.40) and (5.42) of Ref. [H] and starting from 

PT = ^J{pl+P\f + {PI + Plf- (2.29) 



3 NLO results 

We now describe the calculation of the NLO contributions o"o('^), '^ojl^)) ^'^^ a^-^^ of 
Eq. (12. lip in some detail. 

We employ the following tools. We generate the relevant Feynman diagrams using 
the symbolic program package FeynArts [15], carry out the spin and colour sums using 
the program package FormCalc [I6] , and perform the Passarino- Veltman reduction of the 
tensor one-loop integrals [17] using the program package FeynCalc [18]. Subsequently, 
we implement the analytical results in a Fortran program. We evaluate the standard 
scalar one-loop integrals contained in the purely weak corrections using the program 
package Loop Tools [T9], which incorporates the program library FF [20]. The numerical 
integrations are performed using the program package Cuba [21j, which provides several 
different integration routines and is, therefore, also well suited for cross checks. 



3.1 Virtual electroweak corrections to W + j production 

The virtual electroweak corrections of 0{a) to processes (I2.2p - (12.4I) arise from self-energy, 
triangle, box, and counterterm diagrams. They are shown for process ( 12. 2p in Figs. [2HSI 
respectively. 

Evaluating the transition matrix element ^?o(o) from these loop diagrams, we encounter 
both ultraviolet (UV) and IR singularities, which need to be regularised and removed. 
As usual, we use dimensional regularisation, with D = 4 — 2e space-time dimensions and 
't Hooft mass scale fi, to extract the UV singularities as single poles in e. These are 
removed by renormalising the parameters and wave functions of the LO transition matrix 
element 7^^"', which leads to the counterterm contribution (see Fig. [5]), 

n-Wj _ q-Wj^Wj /o i\ 

Owing to the renormalisability of the SM, the UV singularities in '^^^^ cancel, and the 
physical limit e — *• can be reached smoothly. 

The electroweak on-shell renormalisation scheme uses the fine-structure constant a 
defined in the Thomson limit and the physical particle masses as basic parameters. In 
order to avoid the appearance of large logarithms induced by the running of a to the 
electroweak scale Mw in T^^^^, it is useful to replace a by Gi? in the set of basic parameters, 
by substituting 
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where Ar [22] contains those radiative corrections to the muon hfetime which the SM 
introduces on top of those derived in the QED-improved Fermi model. In the electroweak 
on-shell scheme thus modified, we have 



SZ,-^-^ + \ (SZ., 



L 

VM 



w 



Ar) 



(3.3) 



where the renormahsation constants read 



SZ. 



Re 



,2 
w 

2" 

w 



1 c; 

2 s 



w 



-Re 



,2 



Ml 



Ml ^ 



-ReS^,(m^) 



-Ml 
d 



qq 



(3.4) 



Here, S^^^^, S^^, S^^, and E^^ are the transverse parts of the respective electroweak 
gauge-boson self-energies and mixing amplitudes, S^^, E^, and E^^ are the left-handed, 
right-handed, and scalar parts of the quark self-energy, and = 1 — s^. 

The IR singularities can be of soft or collinear type. The loop diagrams involving 
virtual photons interchanged between external lines are plagued by soft IR singularities. 
Owing to the Bloch-Nordsieck theorem |8] , they cancel against similar singularities arising 
from the real emission of soft photons, to be discussed in Section 13.31 The loop diagrams 
involving external quark or antiquark lines that split into virtual photons and quarks 
generate collinear IR singularities. Such singularities also arise from the real emission of 
collinear photons off external quark or antiquark lines, as will be explained in Section 13. 3[ 
Thanks to the Kinoshita-Lee-Nauenberg theorem [9], collinear IR singularities from FSR 
are completely canceled in the sum of real and virtual corrections provided that the final 
state is treated inclusively enough. On the other hand, collinear IR singularities from ISR 
survive and have to be absorbed into the quark and antiquark PDFs. For consistency, the 
splitting functions in the evolution equations of the PDFs then need to be complemented 
by their 0{a) terms. IR singularities also arise from the wave-function renormalisations in 
Eq. (13.41) . We choose to regularise the IR singularities by assigning infinitesimal masses, 
A, niu, and m^, to the photon, the light up-type quarks, and the down-type quarks, 
respectively. This is convenient because the standard scalar one-loop integrals Cq and 
Dq that emerge after the tensor reduction [T7] are well established for this regularisation 
prescription [2l]. Although the purely weak loop corrections are altogether devoid of 
IR singularities, terms logarithmic in rriu and nid are generated by the tensor reduction. 
However, these artificial IR singularities cancel among themselves. 

We emphasise that, in the treatment of both the virtual and real corrections, terms 
depending on A, iriu-, and md are extracted analytically and their cancellation is established 
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manifestly, so that the expressions used for the numerical analysis do not contain these 
IR regulators. 



3.2 Virtual QCD corrections to + 7 production 

The virtual QCD corrections of 0{as) to process (12.51) arise from the self-energy, triangle, 
and box diagrams shown in Fig. [6] and the counterterm contribution. 



(3.5) 



The latter only receives contributions from the gluon-induced wave-function renormalisa- 
tion of the external quark lines. 



(3.6) 



where 



-E?i^(m2) 



o 2 9 
Om 



(3.7) 



Because parity is conserved within QCD, the quark self-energy has just one vector part 
TP^^ = I]^^ = Sg^^. Up to terms that vanish in the limit nig 0, we have 



4:71 



- 7b + ln(47r) -ln^-21n-|+4 



A2 



+ 0{e) 



{3.t 



where Cp = {N^ — l)/{2Nc) =4/3 for Nc = 3 quark colours, 'Je is the Euler-Mascheroni 
constant, and A now represents an infinitesimal gluon mass. 



3.3 Real corrections due to VF + j + 7 production 

The tree-level diagrams for process (12.81) are shown in Fig. [71 They contribute at the same 
time to the electromagnetic bremsstrahlung in process (12. 2p and to the QCD bremsstrah- 
lung in process (12. 5p . which complicates the treatment of the electroweak corrections to 
W+j associated production, as explained in the Introduction. The diagrams contributing 
to the electromagnetic bremsstrahlung in processes (12. 3p and (12. 4p emerge from Fig. [7| by 
crossing the gluon with the u and d quarks, respectively. 

When the cross sections of processes (I2.8p - (l2.10p are integrated over their three- 
particle phase spaces, one encounters IR singularities of both soft and collinear types. 
The former stem from the emission of soft photons and gluons and cancel against simi- 
lar contributions from the virtual corrections owing to the Block-Nordsieck theorem [S], 
as explained in Section 13. 1[ The latter arise when a massless gauge boson is coUinearly 
emitted from an external massless fermion line or when a massless gauge boson splits into 
two collinear massless fermions. Specifically, in process (12.80 . the photon or the gluon can 
be emitted collinearly from the incoming Ui and dj quarks; in process (12. 9p . the photon 
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can be emitted collinearly from the incoming Ui quark or the outgoing dj quark, and 
the gluon can spht into a coUinear djdj quark pair; and in process (12.101) . the photon 
can be emitted colhnearly from the incoming dj quark or the outgoing Ui quark, and the 
gluon can spht into a cohinear UiUi quark pair. As aheady mentioned in Section 13.11 the 
coUinear IR singularities from FSR are canceled by the virtual corrections according to 
the Kinoshita-Lee-Nauenberg theorem [9] if the considered process is treated inclusively 
enough. By contrast, those from ISR survive and have to be absorbed into the PDFs. 

Due to the minimum-p^ cut, the photon and the gluon cannot be soft simultaneously 
because one of them has to balance the transverse momentum of the W boson. By the 
same token, there can only be one coUinear situation at a time. However, soft and coUinear 
singularities do overlap, and care needs to be exercised to avoid double counting. 

For consistency, also the IR singularities in the real corrections need to be regularised 
by the photon and gluon mass A and the light-quark masses ruu and rud introduced in 
Sections 13.11 and 13.21 As already mentioned in Section 13.11 their cancellation is achieved 
analytically, so that the expressions underlying the numerical analysis are free of them. 

As in Ref. |¥|25|26j . we employ the method of phase space slicing [27] to separate the 
soft and coUinear regions of the phase space from the one where the momenta are hard 
and non-coUinear, so that the partonic cross section can be written as 

da^^' = + daZr + daZil (3-9) 

For definiteness, let us assume that parton d in process f l2.12p is the soft or collinearly 
emitted one and that partons a and c are the ones emitting ISR and FSR, respectively. In 
the notation introduced in Section 12.1.11 the soft regions of phase space are then defined 
by A < p2 < AE < (s - w)/{2y/s), the coUinear ones for ISR and FSR by ^9 < A^J < vr 
and ip < Alp ^ tt, respectively, and > AE, and the hard and non-colUnear one by the 
rest. In Sections l3.3.1l and l3.3.2l we explain how to evaluate da^^^"' and do"^j]'^ analytically 
using appropriate approximations. On the other hand, da^Jd can straightforwardly be 
evaluated numerically with high precision [2H]. Since AE is to be measured in units of 
\^/2, we define 6s = 2AE/^/s. The demarcation parameters 6s, Ai), and Atp must 
be chosen judiciously. If the are too small, then the numerical phase-space integration 
performed for da^J^ becomes unstable; if they are too large, the approximations adopted 
for da^l/' and da^J^^ become crude. In practice, one varies 6s, A{}, and Atp to find the 
respective stability regions. For the problem considered here, this is easily achieved. 



3.3.1 Soft singularities 

In the soft phase space regions, 7""^^^cd f^ctorises into 7"'^''^^'= times an eikonal fac- 
tor that depends on p^. Squaring T'^^'^^^d^ performing the spin and colour sums, and 
integrating over pd with the constraint A < p° < AE, one has [231129] 

d^sofT^"^!^, A^) = Cr^"^(A, A^)da"''^^^ (3.10) 
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In the case of soft electromagnetic and QCD bremsstrahlung in process f l2.8p . we then 
obtain 



Colr^ ^^{^^ A-E) — — — {Qu6uu + Qd^dd + SwW + "^QuQd^ud + "^Qu^uW + "^Qd^dw) 



a 



soft V. / 271 



277 



{^uu + ^dd + 25nd) 



(3.11) 



where = 2/3 and Qd = —1/3 are the fractional electric charges of the u and d quarks, 
respectively, and 



5ww 



dW 



In 



A2 



+ ln^ 



^dd — ^uu 
5ud 



2 2 ,1,^ 2 

In \ ^ ^ In 
2 A2 



1 1. 4(AE)2 m^ml ^ 1 ^ !^ ^ 1^2 
2 4 



+ 



TC 



In — + In -, 

s — w s 

1, 4(AE)2, u;m' 



2^"- A^ 



+ Lis 



In 



{w-ty 4 



!L_ + 1 ( ln2 !^ + ln2 ^ 



w — t 



Li2 



W — t 



TT 



-5. 



(3.12) 



Here, Li2(x) = — dt ln(l — tx)/t is the dilogarithm, and terms that vanish for m„ 
iTid = have been omitted. 

Furthermore, we find the soft-photon correction factor for process f l2.9p to be 



^sofT^ '^"'{^■! A-E) — — — ( Qu^uu + Qd^dd + ^WW + '^QuQd^ud + SQu^uVl^ + 2Qrf5(iVK 



in which two terms of Eq. (13. 12^ are modified to be 



(3.13) 



A2 



In^ ^ + In^ ^ 

s [s — wY 



+ Li2 



- I + — , 

u 



~ 1 4(AE)2 wml 1 

Jd^y = -- In In ^ - - 

2 A^ (s — w)^ 4 



In^ — + In 



2 ^"^d 
fs — w) 



Li. 1 



(3.14) 



Finally, the soft-photon correction factor for process (12.10p emerges from the one of 
process (12.91) through the simple replacement 



soft 



''^''"^(A,AE) = SZ7'^^''\\,^E) 



(3.15) 
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3.3.2 Collinear singularities 

As explained in Section [3731 collinear singularities arise from three sources: (1) the emis- 
sion of a photon or gluon from an incoming quark or antiquark; (2) the splitting of an 
incoming gluon into a quark- antiquark pair; and (3) the emission of a photon from an 
outgoing quark or antiquark. The resulting contributions to da^f^ all factorise into the 
respective LO cross sections without radiation and appropriate collinear radiator func- 
tions p^|l30] . In the case of ISR, this also involves a convolution with respect to the 
fraction z of four-momentum that the emitting parton passes on to the one that enters 
the hard interaction. 

Let parton a in process fl2.12p be the emitting quark q and parton d the emitted photon 
or gluon. Then we have [14|30] 



{aQ\ 



asCp] 



2n 



dzR'g'{mg,A^,z) d&Q 



qh^Wc 



Pq-'Zpq 



(3.16) 



where 5s is introduced to exclude a slice of phase space that is both soft and collinear and 

'soft ' ^0 



is already included in aY''"' 



Kim^Ai^^z) = Pg^,{z) 



tqS/s, with To being defined in Eq. fl2.15p . and 

s{A^)^ 2z 



In 



1 + z^ 



(3.17) 



with 



1 + z' 



(3.18) 



being the LO q ^ q splitting function [3TJ. This result readily carries over to the case 
when parton a is an antiquark q. Note that the cm. frame is boosted along the beam 
axis by the collinear emission of the photon or gluon. 

Now, let parton a in process f l2.12p be a gluon that splits into a qq pair, with q being 
outgoing and q entering the residual hard scattering. Then we have ^26j 



m^. Ad) 



2tx 



dzRf{mg,Ad,z) d&l 



Pq^ZPq 



(3.19) 



where Tp = 1/2 and 



Rf{mg,A§,z) 



sil 



Arrig 



+ 2z{l-z), 



(3.20) 



with 



■,{z)=z' + {l-zf 



(3.21) 



being the LO g ^ q splitting function. This result readily carries over to the case when 
parton d is an antiquark q. 
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Finally, let parton c in process (12.121) be the emitting quark q and parton d the emitted 
photon. Then we have [Till30] 



rUg, Alp) 



l-Ss 



dzi?f K,A^,^)dao"'^^^ 



(3.22) 



where 6s = s6s/{s — w) is again to avoid double counting of phase space regions that are 
both soft and collinear, and 



i?f(m„A'(;),z) = P,_ 



2z 



1 + ^2 



(3.23) 



with Pq^q given in Eq. (I3.18p . This result readily carries over to the case when parton c 
is an antiquark q. The integral in Eq. (I3.22p is not a convolution and can easily be carried 
out, yielding 



dzRf{mq,Ai,,z) 



-2\n5s-- I In 



+ 21n5.-3vr + -. (3.24) 



In order to obtain da^j-J'*' for one of the processes (12. 80 - 02. 101) . all possible collinear 
emissions must be taken into account one by one. 

While the collinear IR singularities from FSR cancel upon combination with the virtual 
corrections by the Kinoshita-Lee-Nauenberg theorem [H], those from ISR survive. Since 
their form is universal, they can be factorised and absorbed into the PDFs [32] . Adopting 
the modified minimal-subtraction (MS) factorisation scheme both for the collinear sin- 
gularities of relative orders 0{a) and 0{as), this is achieved by modifying the PDF of 
quark q inside hadron A as 



fq/A{x,M^)^fq/A{x,M' 

- In^ 6s - In 6s + I 
M2 



■■fg/Aix,M 



1 - 



In- 



1 — z)'^m'^ 



- 1 



aQl + agCp 



27r 



2tt 



In (5, + - 

M2 



In 



M2 



Pg^q{z)\n- 



(3.25) 



where M is the factorisation mass scale, which separates the perturbative and non- 
perturbative parts of the hadronic cross section. 



4 Numerical results 

We are now in a position to present our numerical results. We start by specifying our 
choice of input. We adopt the values Gp = 1-6637 x 10"^ GeV-^, Mw = 80.403 GeV, 
Mz = 91.1876 GeV, and = 174.2 GeV recently quoted by the Particle Data Group [55] . 
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take the other Uf = 5 quarks to be massless partons, and assume Mh = 120 GeV, which 
is presently compatible with the direct search limits and the bounds from the electroweak 
precision tests |33j. We take the absolute values of the CKM matrix elements to be |33j 

iKdl = 0.9377, iKsl = 0.2257, = 0.230, 

\Vcs\ = 0.957, iKfel = 41.6 X 10^^, |Kfe| = 4.31 x lO'^. ^ > 

Since we are working at LO in QCD, we employ the one-loop formula for We use 

the LO proton PDF set CTEQ6L1 by the Coordinated Theoretical-Experimental Project 
on QCD (CTEQ) [34j, with Aq^j^ = 165 MeV. In the case of photoproduction, we 
add the photon spectra for elastic [lOj and inelastic [TT1IT2] scattering elaborated in the 
Weizsacker- Williams approximation. In the latter case, we use the more recent set by 
Martin, Roberts, Stirling, and Thorne (MRSTQED04) |I2] as our default, and the set by 
Gliick, Stratmann, and Vogelsang (GSV) pjj to assess the theoretical uncertainty from 
this source. We choose the renormalisation and factorisation scales to be = M = ^m§i^*, 
where m^^^ = ^(p™t)2 + is the minimum transverse mass of the produced W boson 
and is introduced to estimate the theoretical uncertainty. Unless otherwise stated, we 
use the default value ^ = 1. 

We consider the total cross sections of pp — > + X at the Tevatron (run II) with 
= 1.96 TeV and pp ^W^+X at the LHC with VS = U TeV as functions of p^^K By 
numerically differentiating the latter with respect top^*, we also obtain the corresponding 
Pt distributions as dcr/dpT = — da{p'^^) /dp^^l^^t^p^. Owing to the baryon symmetry 
of the initial state, the results for and W~ bosons are identical at the Tevatron, 
and it is sufficient to study one of them. By contrast, iy+-boson production is favoured 
at the LHC because the proton most frequently interacts via a u quark. Therefore, it 
is necessary to study the production of and W~ bosons separately at the LHC. We 
compare the contributions of four different orders: (1) the LO contribution of 0{aas) from 
processes fl2.2l) - fl2.4p . where the system X accompanying the W boson contains a hadron 
jet; (2) the LO contribution of from process (12.51) . where X contains a prompt 

photon; (3) the NLO contribution of O^a'^ag) comprising processes (I2.2p - (I2.5I) at one loop 
as well as processes (I2.8I) - (I2.10I) at tree level, where X contains a hadron jet, a prompt 
photon, or both; and (4) the LO contributions of 0{a^) from processes (12. 6p and (12.71) 
via direct photoproduction and from processes (I2.2p - (l2.4p via resolved photoproduction, 
where X contains a hadron jet and, in the case of elastic photoproduction, also the 
scattered proton or antiproton. Since we consider inclusive one-particle production, we 
do not use any information on the composition of X, i.e. we include all possibilities. In the 
following, we regard the sum of contributions (1) and (2) as LO and sum of contributions 
(l)-(4) as NLO unless the perturbative orders are explicitly specified in terms of coupling 
constants. We thus define the correction factor K to be the NLO to LO ratio with this 
understanding. 

Let us now discuss the numerical results and their phenomenological implications in 
detail. Specifically, Figs. [HI Elta), and [TUlfa) refer to the Tevatron, while Figs IWb). [TUlfb). 
[TTl [T2l and [13] refer to the LHC. In Fig. [8](a) the NLO result for the total cross section 
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as a function of p™*^ is compared with the LO contributions of 0{aas) and 0{a^) as 
well as with the photoproduction contribution of 0{a^). The 0{aas) and C(a^) results 
exhibit very similar line shapes, but the normalisation of the latter is suppressed by a 
factor of about 500. This may be qualitatively understood from the partonic cross section 
formulae in Eq. fl2.2ip and by noticing that the 0{aas) contributions from the Compton- 
like processes ( 12.31) and (12. 4p . which have no counterparts in O(a^), are significantly 
enhanced by the gluon PDF. As a consequence, the LO result is almost entirely exhausted 
by the 0{aas) contribution. 

The inclusion of the NLO correction leads to a moderate reduction in cross section, 
which increases in magnitude with p™*, reaching about —4% for p™^ = 200 GeV, as may 
be seen from Fig. M^a), where the K factor is depicted. 

In Fig. [8]^a), also the photoproduction contribution is shown. As explained above, we 
have to distinguish between elastic and inelastic scattering off the proton on the one hand, 
and between direct and resolved photons on the other hand, so that, altogether, we have 
four different contributions, which all formally contribute at 0{a^). The resolved-photon 
contributions turn out to be small against the direct-photon ones and are, therefore, 
not included in Fig. M^a). As for the combined direct-photoproduction contribution, we 
observe from Fig. [8]^a), that, except for small values of p™*, it overshoots the O(a^) contri- 
bution, although it is formally suppressed by one power of a\ Detailed inspection reveals 
that this unexpected enhancement can be traced to the direct-photoproduction diagram 
involving the triple-gauge-boson coupling and the space-like W^-boson exchange, which 
significantly contributes at large values of ^/s. In fact, for a fixed value of p™*, the total 
cross sections of processes (12.61) and (12. 7p have an asymptotic large-s behaviour propor- 
tional to l/(m™*)^, while those of processes (12. 20 - 02. 50 behave as Ins/s. Consequently, 
photoproduction appreciably contributes to the K factor, as is apparent from Fig. M^a), 
which also shows the photoproduction to LO ratios for elastic and inelastic scattering. 
The freedom in the choice of the inelastic photon content of the proton is likely to be the 
largest source of theoretical uncertainty in the photoproduction cross section. In order to 
get an idea of this uncertainty, we display in Fig. M^a) also the inelastic-photoproduction 
to LO ratio evaluated with the GSV photon spectrum for inelastic scattering. The result 
is roughly a factor of two smaller than our default prediction based on the MRSTQED04 
spectrum. 

In Fig. [TOl(a). we examine the theoretical uncertainties in the 0{aas), C(a^), NLO, 
and photoproduction results due to the freedom in setting the renormalisation and factori- 
sation scales by exhibiting their ^ dependencies relative to their default values at ^ = 1. 
The ^ dependencies of the 0{a'^) and direct-photoproduction results stem solely from 
the factorisation scale M and are rather feeble, while those of the 0{aas) and resolved- 
photoproduction results are also linked to the renormalisation scale /i of as{n) and are 
more pronounced, but still not dramatic. The scale variation of the LO result amounts to 
less than ±15% for 1/2 < ^ < 2. It is only slightly reduced by the inclusion of the NLO 
correction. This is expected because the NLO result is still linear in ^^(/i), so that the n 
dependence of ^^(/i) is not compensated yet. 

In Fig. Mjo), the analysis of Fig. [8](a) is repeated for the pt distribution. We observe 
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that the hne shapes and relative normahsations of the various distributions are very similar 
to those in Figs. M^a) and the same comments apply. 

Turning to the LHC, we can essentially repeat the above discussion for the Tevatron, 
except that we have to take into account the difference between and W~ boson 
production. Thus, Fig. M has two LHC counterparts, Figs. [H] and [121 for the and 
W~ bosons, respectively. To illustrate this difference more explicitly, we show in Fig. [13] 
the to W~ ratios of the respective results from Figs. [11] and [12] For simplicity, 
Figs. [9](b) and[TO](b). the LHC counterparts of Figs. [9](a) and [TCT a). refer to the averages 
of the results for and W~ bosons. In the following, we only focus on those features 
which are specific for the LHC. From Figs. [11] and [121 we observe that the gaps between 
the 0{aas) and results are increased by about a factor of two, to reach three 

orders of magnitude. This is mainly because the Compton-like processes (12.31) and (12.41) 
benefit from the extended dominance of the gluon PDF at small values of x. Furthermore, 
the photoproduction contributions now significantly exceed the C(a^) ones throughout 
the entire p™* and pt ranges. From Fig. [T31 we see that the to W~ ratios take 
values in excess of unity, as expected, and strongly increase with increasing values of 
p™*. Comparing Figs. [9](a) and (b), we find that the K factors are significantly amplified 
as one passes from the Tevatron to the LHC. This is due to the fact that the Sudakov 
logarithms, which originate from triangle and box diagrams, become quite sizeable at the 
large values of and pt that can be reached at the LHC. This issue was already dwelled 
on in Ref . [7] , to which we refer the interested reader. Finally, comparing Fig. [TOT a) and 
(b), we conclude that the ^ dependence is generally somewhat smaller at the LHC. 

5 Conclusions 

We studied the effect of electroweak radiative corrections at first order on the cross section 
of the inclusive hadroproduction of single W bosons with finite values of px, putting special 
emphasis on the notion of infrared-save observables with a democratic treatment of hadron 
jets initiated by (anti) quarks and gluons. This is indispensable b matter of 

principle, a coUinear gluon-photon system cannot be distinguished from a single gluon 
with the same momentum, so that a minimum-transverse-momentum cut on the gluon is 
an inadequate tool to prevent a soft-gluon singularity. This led us to include the 0{as) 
correction to W + •y production along with the 0{a) correction to W + j production, 
both contributing at absolute order 0{a'^as). We also considered the contribution from 
events where one of the colliding hadrons interacts via a real photon, which is of absolute 
order 0{a^). The hadron can then either stay intact (elastic scattering) or be destroyed 
(inelastic scattering), and the photon can participate in the hard scattering directly (direct 
photoproduction) or via its quark and gluon content (resolved photoproduction), so that 
four combinations are possible. 

We extracted the UV singularities using dimensional regularisation and removed them 
by renormalisation in the on-shell scheme. We regularised the soft and coUinear IR sin- 
gularities by means of infinitesimal photon, gluon, and quark masses. A, m„, and m^. 
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respectively. We used the phase-space shcing method, with cuts Sg, Ai?, and A'^ on the 
scaled photon and gluon energies and on the separation angles in the initial and final 
states, respectively, to isolate the soft and collinear singularities within the corrections 
from real particle radiation. We achieved the cancellation of A, m„, and nid analytically 
and ensured that the numerical results are insensitive to variations of Ss, At?, and A'^ 
within wide ranges about their selected values. 

Wc presented theoretical predictions for the total cross sections with a minimum-p2- 
cut and for the px distributions to be measured in pp collisions with = 1.96 TeV 
at run II at the Tevatron and in pp collisions with = 14 TeV at the LHC, and 
estimated the theoretical uncertainties from the scale setting ambiguities. We found that 
considerably less than 1% of all VF + X events contain a prompt photon. The corrections 
considered turned out to be negative and to increase in magnitude with the value of Pt- 
While the reduction is moderate at the Tevatron, reaching about —4% at Pt = 200 GcV, 
it can be quite sizeable at the LHC, of order —30% at pt = 2 TeV, which is due to the 
well-known enhancement by Sudakov logarithms. It is an interesting new finding that 
the photoproduction contribution is considerably larger than expected from the formal 
order of couplings. In fact, it compensates an appreciable part of the reduction due to 
the 0{a^as) correction. 
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(b) 

Figure 1: Tree-level diagrams of (a) process fl2.2p and (b) process fl2.5l) . The tree-level 
diagrams of processes (12. 3p . fl2.4p . fl2.6p . and fl2.7p emerge through crossing. 
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Figure 2: 0{a) self-energy diagrams of process (12.21) . 
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Figure 6: 0{as) self-energy, triangle, and box diagrams of process fl2.5p . 




Figure 7: Tree- level diagrams of process (12. 8p . The tree-level diagrams of processes 
and fl2.10p emerge through crossing. 
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Figure 8: (a) Total cross section as a function of p™* and (b) pt distribution of pp — > 
+X for = 1.96 TeV (Tevatron run II). The NLO results are compared with those 
of orders 0{aas), C(a^), and 0{a^) via photoproduction. 
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Figure 9: NLO corrections [K — 1), with and without the photoproduction contributions, 
to the total cross sections of (a) pp + X for \fS = 1.96 TeV (Tevatron run II) and 

of (b) pp^W^+Xhi^= 14 TeV (LHC) as functions p^^\ For comparison, also the 
contributions due to elastic and inelastic photoproduction normahsed to the LO results 
are shown. In the latter case, the evaluation is also performed with the GSV PDFs. 
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Figure 10: Total cross sections of (a) pp W++X for = 1.96 TeV and p^^^ = 20 GeV 
(Tevatron run II) and of (b) pp + X ioi VS = 14 TeV and p"^^ = 200 GeV (LHC) 
as functions of ^ normalised to their default values for ^ = 1. The NLO results are 
compared with those of orders 0{aas), C'(a^), and O(a^) via photoproduction. 
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Figure 11: (a) Total cross section as a function of p^^ and (b) pt distribution of pp — > 
+ X for = 14 TeV (LHC). The NLO results are compared with those of orders 
0{aas), O^a"^), and 0{a^) via photoproduction. 
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Figure 12: (a) Total cross section as a function of p^^* and (b) px distribution of pp — > 
W~ + X for = 14 TeV (LHC). The NLO results are compared with those of orders 
0{aas), Oi^a^), and Oi^ct') via photoproduction. 
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Figure 13: Ratios of the respective results for and W bosons shown in Figs. [TlT a) 



27 



